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Abstract 

We introduce a "pre-Bethe-Ansatz" system of equations for three dimensional vertex mod- 
els. We bring to the light various algebraic curves of high genus and discuss some situations 
where these curves simplify. As a result we describe remarkable subvarieties of the space of 
parameters. 
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1 Introduction 



The purpose of this paper is to sketch how ideas introduced in the study of the sixteen vertex 
model in can be generalized to higher lattice dimensions. We think that the ideas developed 
here are relevant tools for the analysis of lattice models in three or more dimensions, a widely 
unexplored area. 

In this paper we first recall the basic results obtained for the two-dimensional case, with a special 
emphasis on the symmetries they uncover. We then introduce the simplest three-dimensional 
generalization of these results. In order to be more concrete, we describe a specific model which 
naturally generalizes the Baxter model. Finally, we show how this general construction points to 
a number of algebraic varieties of interest. 



2 Preliminaries to the Bethe Ansatz 

One of the keys to the Bethe Ansatz for two-dimensional vertex models is the existence (see for 
instance equations (B.10), (B.lla) in Q) of vectors which are pure tensor products (of the form 
wgiui), and which the i?-matrix maps onto another pure tensor product v' <8>w' ||. The "pre-Bethe 
Ansatz" equation is therefore: 

R (v ® w) = fJ, v' (g> w', (1) 

A group G Bethe — sl 2 x sl 2 x sl 2 x sl 2 naturally acts on this equation, extending the weak-graph 
transformations |l| . With the following notations 




the elimination of w and w' from (|l]) yields the biquadratic relation: 

aoo + aio P + oe i p 1 + a 20 p 2 + a 02 p' 2 + 2an pp' + a 2 \ p 2 p' + a i2 pp' 2 + a 22 p 2 p' 2 = 0. (3) 

The ay's are quadratic polynomials in the homogeneous parameters ai, • • • , dj of the i?-matrix. 

Equation (^) (and the similar one for q and q') is a quadratic Frobenius relation ^|, |^| on 
theta functions. Let us recall that the Frobenius relations are intertwining relations of the product 
of two theta functions Q. The "intertwining" matrix R is expressible in terms of ratio and 
product of theta functions depending on the difference of the arguments (spectral parameters) of 
the theta functions appearing in the parametrization of p, p' (resp. q, q') Qj. This intertwining 
functional relations on theta functions enable to represent a Zamolodchikov algebra (see Q) which 
is "almost" [] a sufficient condition for the Yang-Baxter equations to be satisfied: 

Mb) Mh) = - 9 ^ Ak ^ MOi) (4) 

kl 

In [Q, it has been shown that equation (^) defines an elliptic curve in the p,p' plane. The 
discriminant of (g) seen as a quadratic polynomial in p is: 

D = ap' 4 + 4f3p' 3 + 6jp' 2 +A/3'p' + a' . (5) 

The transformations of D under the sl 2 group acting on v' (homographic transformations of p') 
yield two fundamental invariants g 2 and 53 |9|, [l(| and a modular invariant J = g\j (g 2 — 27<?|): 

g 2 = aal -4(3/3' + 3 7 2 (6) 
g 3 = aja' + 2(3jf3' - a/3' 2 - a' (3 2 - 7 3 (7) 

1 The very existence of a Zamolodchikov algebra implies the Yang-Baxter equations to be satisfied provided the 
linear independence of the expressions Ai(jh)Aj (02)-Ak {63)- This independence condition is apparently not satisfied 
for theta functions of several parameters Jm. 
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Exchanging the role of p andp' in equation (||) leads to the same 172, <?3 and J. More remarkably, 
the same operations applied to the equation in q and q' similar to (|^), lead to the same gi and 
53. This is a quite non-trivial result. One can therefore associate the same Weierstrass's canonical 
form (elliptic curve), y 2 = Ax 3 — g2X — 33, to equation (|3|), the similar equation in q and q' and 
also to the elliptic curve, intersection of quadrics in CP15, which corresponds to the orbits of T 
(see Figure 1 in [Q). 

The algebraic variety defined by J = 00, i.e. g\ — 27 g 2 = 0, is a remarkable variety in the 
parameter space CP15: the elliptic parametrization of ^) becomes a rational one. Recalling the 
well-known example of the symmetric eight-vertex Baxter model Jll], [l2|], this reduction to a 
rational parametrization corresponds to two remarkable situations in the parameter space: the 
disorder varieties, for which some dimensional reduction occurs fl3|| , and the critical varieties fli] , 
[l]. Quite surprisingly, these two varieties of a priori different nature appear here on the same 
footing. It would be pointless to write down the condition J = 00, since it is a sum of several 
millions of monomials of degree 24 in the sixteen homogeneous parameters of the sixteen- vertex 
model. 

Noticeably, the Frobenius relations do exist for theta functions of several variables (g variables) 
leading to some associated Zamolodchikov algebra Q. However, they do not lead to Yang-Baxter 
equations J8J . Clearly, there is room for models satisfying a non- trivial functional relation like 
equation (|l|) (which is a key for the explicit construction of the Bethe Ansatz) but not the Yang- 
Baxter equations. We may for instance consider the generalization of the symmetric eight- vertex 
model constructed in Q for which there exists a parametrization in terms of theta functions of g 
variables: this model can be seen as a coupling of g replicas of the Baxter model, the arrows of the 
vertex taking values in Ti?, and the row-to-row transfer matrix (with periodic boundary conditions) 
being a 2 gN x 2 gN matrix. If the Yang-Baxter equations were satisfied (which is apparently not 
the case [||), they would yield the existence of a ((/-dimensional) family of commuting transfer 
matrices. This is certainly too demanding. The existence of transfer matrices commuting only 
in a subspace of the 2 gN vector space, corresponds to interesting models and may be sufficient 
to calculate their partition function (largest eigenvalue of the transfer matrix). One must recall 
the example of the disorder solutions [l^] , for which a commutation of the (diagonal) transfer 
matrices is satisfied on a one-dimensional vector space, which is actually sufficient to calculate the 
partition function (see Jl5| and page 247 of @). 

The exchange of p and p' (resp. q and q'), as well as the symmetries described in Jlql , are 
realized as shifts of the spectral parameter 9 of the elliptic curve (see for instance |]l6| ). This is a 
key point for the explicit construction of the Bethe Ansatz, since it enables to build a vector 

!*(*)) =(g)e(0 + 277n), (8) 

n 

which transforms very simply by the row-to-row transfer matrix: 

|*(0)>-|*(0 + 2»,)>. (9) 

When acting on ^, the translation by a lattice spacing amounts to the shift of the spectral 
parameter 9. Taking into account the translation invariance of the lattice with periodic boundary 
conditions, one can "Fourier transform" | and get an eigenvector \^{9)) of the transfer matrix 
(see p. 148 of |l6| ) . Notice that this construction works for theta functions of g variables using the 
corresponding Frobenius relations One therefore has a commutation of transfer matrices on 
this (one dimensional !) vector space. This is just the first step for the construction of the Bethe 
Ansatz. The second step amounts to build eigenvectors made up of product vectors (similar to 
(||)) with a given number of "deviating" factors (see p. 148-150 of |ll| and ||). 
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3 Towards three-dimensional Bethe Ansatz 



A similar analysis can be performed in higher dimensions. In this section, we show how we can 
generalize these results in dimension three. In a specific example, we find three curves linked to a 
matrix R living on a higher dimensional variety | p7[ . 

We denote w(i, j, k,l,m, n) the Boltzmann weight of a given three-dimensional vertex. We shall 
only consider the simplest case where each of the spins i, j, k, I, m and n can take only two values. 
The vertex weights may be arranged in an 8 x 8 matrix of entries: 

R Znn =W ( i '3,k,l,m,ri). ( 10 ) 

The natural generalization of the "pre-Bethe Ansatz" equation ([!]) is: 

R (u <g> v ® w) = /J, ■ u' ® v' (g> w'. (11) 

Let us introduce the notation: 




In the following subsections, we shall recall the symmetries of a three dimensional vertex model, 
as described in jl^j . 



3.1 The group of inversions T^r, 

As in |is[| , we first introduce the involution I changing R to its matrix inverse (we let appear an 
overall factor A since the entry of R are taken projectively). 

E m%£<» ( 13 ) 



Multiplying both sides of (11) by IR, we get an equation of the same form as ( |ll| ) with u and u', 
v and v 1 and w and w' exchanged and R replaced by IR. 

In [jTsj , we also introduced three partial transposition ti, t 2 and ^3. t\ is defined by: 

The definitions of t2 and £3 are similar. 

The four involutions I and ti (i = 1, 2, 3) generate an infinite discrete group T^d |l8|]. The 
so-called inversion relations of the statistical mechanics model can be simply expressed with these 
building blocks. They are: 

/, J = t 1 It 2 t 3 , K = t 2 It 3 t 1 , L=t 3 Iht 2 . (15) 

Considering the parameter space as a projective space (the entries of the i?-matrix are homogeneous 
parameters), the elements of the group T^d have a non-linear representation in terms of birational 
transformations. This group of symmetry of the parameter space of the model is very large. The 
number of elements of length I grows exponentially with I. It is actually a hyperbolic Coxeter 
group |l9| |. The symmetry group of the Yang-Baxter equations in two dimensions is a mere affine 
Coxeter group |§ |(| (Of. 

The group T^rj has been shown in fl§| to enter the description of the group of automorphisms 
of the tetrahedron equations (generalization of the Yang-Baxter equations in three dimensions). 
We shall use this symmetry group beyond integrability, that is to say for models which do not have 
to verify the tetrahedron equations. 



3 



3.2 Weak- graph duality for 3D models: the gauge group G 

A "gauge" group G = SI2 X SI2 X SI2 acts linearly on the matrix R by similarity transformations 
(the weak-graph transformations, see 0] for details). If g = g\ x 32 x 53, we define: 

g(R) = gi 92 53 ■ R-g^g^gz 1 - (16) 

Each of the gi 's acts on the corresponding vector space and gi for example is a short hand notation 
for gi ® I® I. The action of G and Y^d do not commute. However, G and / do commute, and the 
commutation relation between the t^s and G gives a rather simple semi-direct product structure 
to the combined group: 

tig = g ll h, (17) 

with: 

g U = x 32 x g 3 , (18) 

and similar relations for t2 and £3. In particular, Y^d sends orbits of G onto orbits of G. The 
compatibility of these two groups is described in [fij in a two-dimensional case, the sixteen- vertex 
model. 

The effect of such a transformation on the pre-Bethe-Ansatz equation ([ll]) is simple: g\ act 
naturally on u and it', 32 on u and t/ and 33 on w and u/. 



4 A three-dimensional model 

The most general vertex models on a cubic lattice has a large number of parameters (sixty- four). 
We therefore impose some relations on the Boltzmann weights of the three-dimensional vertex. 
We require that these relations are invariant under the inverse / p2 , |23| and the three partial 
transpositions t\, £2 and t% (equation (|l4|)). They will thus be invariant under the group Y^d- 
We are particularly interested in generalizations of the Baxter symmetric eight- vertex model, and 
define here a specific three-dimensional model, denoted in the sequel B^d- It is possible to "project" 
down a three-dimensional model onto a bidimensional one by just taking the trace of the matrix 
R on one of the spaces 1,2, or 3: take for example space 3. 



n ki 



ET>i,j,ot. 



(19) 



<>3 



The constraints verified by B%d are such that the three possible projections are symmetric Baxter 
models. 

We define &zd by imposing the following restrictions on the entries Eg: 



313233 

3lh33 



n— U>— 12)— »3 
-31,-32,-33 

if hhhhhfa 



-1 



(20) 
(21) 



These constraints imply that the 8x8 matrix R is the direct product of two times the same 4x4 
submatrix p7| . It is further possible to impose that this 4x4 matrix is symmetric, since such a 
symmetry is preserved by the partial transpositions ii, t2, £3 |l^, [l7| ], that is: 



0*1*2*3 _ nJU2j3 

J1J2J3 jl*2*3 



(22) 



We shall use the following notations for the entries of this 4x4 submatrix: 



/ 


ft 




<k 


d 3 \ 




di 


61 


C3 


C2 






C3 


&2 


Cl 


V 


4 


C2 


Cl 


&3 / 



(23) 
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The four rows and columns of this matrix correspond to the states (+, +, +), (+, — , — ), (— , +, — ) 
and (—,—,+) of the triplets (ii, i 2 ,i 3 ) or (ji, j 2 , j 3 )- The matrix R can be completed by spin 
reversal, according to (p0|). Note that t\ (resp. t 2 , ^-3) simply exchanges C2 with d 2 and C3 with d 3 
(resp. circular permutations). / acts as the inversion of this 4x4 matrix. 

It is quite remarkable that there exist four quantities which are covariant by all the four gen- 
erating involutions /, t\, t 2 , £3, and therefore the whole group T^d. Let us introduce: 

Ps = ab 3 + hh - c\ - d 2 3 , q 3 = cidi ~ c 2 d 2 , (24) 

and the polynomials obtained by permutations of 1, 2 and 3. They form a five dimensional space of 
polynomials. Any ratio of these polynomials is invariant under all the four generating involutions 
/, fi, t 2 , £3. CPg is thus foliated by five dimensional algebraic varieties invariant under the whole 



group T 3B . We can also express it by saying that the polynomials (24) define a map from the 
parameter space CPg to CP4 invariant under T 3B) . 

If we consider a subgroup I^ generated by only two involutions, say / and L ( |l5|) or equivalently 
/ and £3 , one gets three more independent covariant polynomials leading to algebraic surfaces (see 
Fig. 1). They read: 

r 3 = ab 3 - bib 2 - cl + dl, 
S3 = (a + b 3 )c3-did2-ciC2, 

t 3 = {bi + b 2 )d 3 - dic 2 - ad 2 . (25) 

From the projection ( |l9| ) we get a Baxter model. If we denote by a Bl b Bl c Bl the non zero 
entries of the i?-matrix of this model, we have: 

a B =a + b 3 , b B = b 1 + b 2 , c B = 2c 3 , d B = 2d 3 . (26) 



5 Study of the three-dimensional "pre-Bethe" equations. 
5.1 A first attempt. 

In the study of the equation (fy]), we can start by eliminating the variables p and p' . We obtain 
the following system of five equations on the remaining variables q, q' , r and r'\ 

= c 3 d 3 - b 2 b 3 q 2 - c\ r 2 + (^3 - b 2 c x + c 2 c 3 - b 3 ci) rq + c 2 d 2 q 2 r 2 - ab x q' 2 
+ (ab 3 - C3 + hb 2 - d\) qq' + (aci - c 2 c 3 + 61C1 - d 2 d 3 ) rq' 
+ (b 2 di - c 3 d 2 - c 2 d 3 + b 3 di) rq 2 q' + 2{c\di - c 2 d 2 ) r 2 qq' 

- (adi + bidi - c 2 d 3 - c 3 d 2 ) rqq' 2 + c 3 d 3 q 2 q 2 + c 2 d 2 r 2 q' 2 - d\ q 2 r 2 q' 2 , (27) 

= {cic 3 - b 2 c 2 ) q + (b 3 c 3 - cic 2 ) r + (ad 2 - b 2 d 3 ) q 2 r + (b 3 d 2 - cxd 3 ) qr 2 
+ ac 2 q' - ac 3 r' - cid 3 q 2 q' - b 3 d 2 r 2 q' 
+ b 2 d 3 q 2 r' + cid 2 r 2 r' + d\d 3 q 2 r 2 q' - did 2 q 2 r r' 

+ (ad 3 + c 2 di — b 3 d 3 - cid 2 ) qrq - (ad 2 + c 3 di - b 2 d 2 - C\d 3 ) qrr , (28) 

= c 3 d 2 - b 2 ci (q 2 + r 2 ) + {c 2 -b 2 + d 2 - c\) rq + (ab 2 - d 2 ) rq' + (aci - ^3) qq' 

+ (bib 2 - C3) qr' - abi q'r' + c 3 d 2 q 2 r 2 + (b 2 di - c 3 d 2 ) r 2 qq' + (cidi - c 3 c?3) rq 2 q' 

+ (cidi - c 2 c?2) r 2 qr' - d\ q 2 r 2 q'r' + (c 2 c?3 + C3^2 - a^i - 61^1) rqr'q' 

+ (61C1 - c 2 c 3 ) rr' + c 2 d2 r 2 qr + (b 2 di - c 3 d 2 ) r'rq 2 + c 3 d 3 q 2 q'r', (29) 
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= c 3 di - b 2 c 2 q 2 - C1C3 r 2 + (did 2 - C\C 2 - b 2 c 3 + &1C3) rq + bid 2 q 2 r 2 
+ (ac 2 - did 3 ) qq' + (ac 3 - d\d 2 ) rq' - ac 3 q'r' 

+ (c 3 d\ - b\d 2 ) r 2 qq' + (c 2 di - bid 3 ) rq 2 q' + c x d 2 r 2 q'r' + b 2 d 3 q 2 q'r' 

- (ad 2 + c 3 di — cid 3 - b 2 d 2 ) rqr'q' - d\d 2 q 2 r 2 q'r' , (30) 

= (c 3 d 2 - b 2 di) q + (c 3 d 3 - cidi) r - (ab 2 - d 2 ) q 2 r - (aci - d 2 d 3 ) qr 2 + ad\ q' 
+ (a 2 + d\ - dl - d\) rqq' - d 2 d 3 (q 2 + r 2 )q' + ad x q 2 r 2 q' + {b x b 2 - c 2 ) qq'r' 
+ (^d - c 2 c 3 ) rr'q' - {ad 1 + b x d x - c 2 d 3 - c 3 d 2 ) rqr'q' 2 + c 2 d 2 r 2 q' 2 r' - ab l q' 2 r' 
+ {c\di - c 2 d 2 ) r 2 qr'q + (b 2 di - c 3 d 2 ) rq 2 r'q - d\ q 2 r 2 q 2 r + c 3 d 3 q 2 q' 2 r' . (31) 

Two similar system of equations can be obtained by the elimination of the pair of variables q and q' 
or r and r' . The equations are of degree two in each of the variables, that is a overall maximum 
degree of eight. In fact only one is of degree seven, three are of degree six and one is of degree five. 

The only apparent property of this system is the invariance by changing the sign of each of the 
variables q, q' , r and r', which is linked to the zeroes of the R matrix for ijklmn = — 1. 

The spin reversal symmetry of the R matrix and the change of R into its inverse IR have no 
visible consequences. This is due to the particular choice made in the elimination of p and p' . The 
equations we have written are just five out of a system of thirty-six equations with a number of 
relations among them. On this full system, the symmetries should be more manifest. 



5.2 Necessary conditions for B373. 

Since this direct attempt to find a full solution to ( pT| ) leads to such a confuse result, we shall 
a study some necessary conditions. We replace the tensor product u ® v and v! ® v' by general 
vectors in the tensor product of space 1 and 2 U and U'. (011) becomes: 



R (C7® w) = nU'®w'. (32) 

This can be written: 

Ri R 2 \ { U\ ( U' \ . . 

M jtj, , (33) 



R 3 R 4 J \rU J r \r'U 

with Ri , • ■ ■ , i? 4 the 4x4 blocs of R. 

Eliminating U and U' gives the necessary condition 



det(Ri r' + R 2 rr' - R 3 -R 4 r)= 0. (34) 



For the model introduced in section (^) (eqs. (|20j), (|2l|), fl22Q), equation (]34|) reads in terms of r 
and r'\ 

A 3 ■ (r 4 r' 4 + 1) + B 3 - (r 4 r' 2 + r 2 r' 4 + r 2 + r' 2 ) + C 3 ■ (r 4 + r' 4 ) 

+D 3 -(r 3 r' 3 +rr') + E 3 -(r 3 r' + rr' 3 ) + F 3 ■ r 2 r' 2 = 0. (35) 

Here A 3 ,...,F 3 are polynomial expressions of degree four in the homogeneous entries of the R- 
matrix ( p3| ) (a, • • • , d 3 ). The simplest expressions are: 

A 3 = (cidx - c 2 d 2 ) 2 = ql, 

C 3 = (b 1 b 2 -c 2 )(ab 3 -d 2 3 ) = ^(p 3 + r 3 )( P3 -r 3 ). 

In fact, all the coefficients A 3l ■ ■ ■ , F 3 can be expressed as quadratic expressions in the polynomials 
invariant by the subgroup Y 2 of T 3 d listed in (24) and (^5|). This shows that equation ( |35| ) is 
invariant by this infinite group. Moreover, they verify the relation: 

4A 3 ■ (F a - 2A 3 + 2C 3 + 2E 3 ) = (D 3 + 2B 3 f . (36) 
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Relation (|3£|) is actually the condition for relation ( |35| ) to become, when r — r' , the square of 

(c 2 d 2 - cidi){r 4 + 1) + r 2 (ab 2 + b x b 3 -c\-d\- ah - b 2 b 3 + c\ + d\). (37) 

We recognize q 3 in the coefficient of r 4 + 1 and p 2 —p% (|24|) in the coefficient of r 2 . These coefficients 
are thus r 3 rj-covariant polynomials. 

Of course two similar eliminations can be performed on (|ff] ) yielding constraints between p and 
p' (resp. q and q'). 

To take into account the symmetry of ( |35| ) by the exchanges r <-> r' , r <-> f/r and r' <-> 1/r', 
one may introduce the variables X = rr' + 1/rr' and Y = r / V + r'/r. Equation (|35|) then becomes 
a conic: 

A 3 X 2 + B 3 XY + C 3 Y 2 + D 3 X + E 3 Y + F 3 = 0, (38) 

with F 3 = F 3 - 2A 3 - 2C 3 . 

An invariant Iproj is naturally associated to the conic (|3S|) , it is the determinant of the 3x3 
matrix M [§: 

/ A 3 S 3 /2 D 3 /2\ 
M = B3/2 C 3 E 3 J2 (39) 
VA3/2 ^ 3 /2 F 3 / 

The value of this invariant is (taking into account the relations between the entries of M): 2p ro j = 
-Z 2 /4^ 3 , with 

1 = 2 B% + D 3 B 3 - 2 E 3 A 3 - 8 A 3 C 3 . (40) 

2proj = is a projectively invariant condition for the conic, meaning that it is the union of two 
lines. Remark that this does not imply the existence of a rational uniformization of (|35|). A similar 
phenomenon happens in the Baxter model, for which (|J) becomes linear in X and Y, but this does 
not provide the elliptic parametrization. 

In order to obtain a parametrization of (^), we look at it as a polynomial of degree four in r, 
the coefficients being polynomials in r': 

a{r')r 4 + 4(3{r')r 3 + 6 7 (r> 2 + 4/3' (r> + a'(r') = (41) 

Its discriminant reads || |f0| |: 

A(r') = 52 (r') 3 - 27. 93 (r') 2 , (42) 

with: 

g 2 (r') = a(r')a'(r')-4/3(r')/3'(r') + 37(r') 2 , (43) 

and: 

g 3 (r>) = a{r') 1 {r')a , {r')+ 2 (3(r') 7 (r'){3'(r') 

-a(r')/3'(r') 2 - /3(r')V(r') - 7 (r') 3 . (44) 

For a general vertex model in three dimensions, equation (|34|) leads to an equation like ( pf| ) 
where a(r'), 0(r'), j{r'), j3'(r'), a'(r') are polynomials of degree four in r' . The polynomials 
92(1"'), g 3 (r'), and A(r') are polynomials of degree 8, 12 and 24 respectively. For the general 
vertex model the hyperelliptic curve y 2 = A(r') is a genus eleven curve. 

However, one verifies easily that, for the model B 3 d, A(r') is a polynomial of degree twelve 
in r' 2 . Moreover the polynomial A(r')/r' 12 is symmetric under the inversion r' <-* f/r'. Hence, 
introducing the variable s' = r' 2 + r'~ 2 , A(r')/r' 12 becomes a degree six polynomial in s', denoted 

If Pq{s') were a generic polynomial of degree six, the hyperelliptic curve y 2 = Pe(s') would be a 
genus two curve, meaning that, as far as parametrization is concerned, one is obliged to deal with 
theta functions of two variables (the Jacobian associated to the genus two curve (6|) or automorphic 
functions (see and page 455 of pj|). One can envisage handy parametrization when the 
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hyperelliptic curves degenerate into elliptic ones, that is, when two roots of Peis') coincide, or 
equivalently when the discriminant of -Pe( s ') vanishes^. 

It is important to note that the model of section (|4|), corresponds to such a situation where 
Peis') can be written as: 

P 6 (s') = (s'-s ) 2 -P i (s') (45) 

where Pa(s') is a polynomial of the fourth degree in s', containing 289 monomials of degree eight 
in the coefficients A.3 , • • • , E 3 . Noticeably, sq is a quite simple expression: 

D 3 + 2B 3 

80 = — YIT (46) 

which reads in terms of the entries a, • ■ • , d 3 of R: 

(cf -c 2 + d 2 - d 2 ) - (a - 63) (61 - h) 

■so = -T 1 ~, (47) 

ci di - c 2 d 2 

Expression fl47| ) is invariant under the group T 3 o (see equations In equation ( |35|), index 3, 

and the equations similar to ( |35| ) relating p and p', or q and q' , lead to equations like (|47|), where 

1.2 and 3 are permuted. It would be interesting to look for conditions on the entries of the R such 
that these three elliptic curves identify. 

Clearly, a particular variety plays a special role: the subvariety in the space of models where 
the three elliptic curves reduce to rational ones. This algebraic variety is a good candidate for 
being a set of critical points (or disorder points) for B 3 d though it is only a codimension-three 
subvariety of the codimcnsion-one critical manifold we are looking for. 

To sum up, the three-dimensional vertex model B 3 rj yields a generalization of the intertwining 
quadratic Frobenius relations in the form of an intertwining of three different elliptic curves by an 
i?-matrix living on an algebraic variety of dimension five given by the intersection of five quadrics 
(equations (|24|)). 

5.3 Further analysis 

The polynomial Pa(s') appearing in equation ( pl5| ) is worth analyzing. We see that 172 and A = 
g\ — 27<7§ factorize: 

g 2 = Al-g ( i ) -g ( i\ g 3 = A$ ■ , A = A% ■ A 1 A 2 A 3 A 4 Ai (48) 



with: 

gi l) = 3 El A z + 8 C3E3A3 + 16 AlC z - D 2 3 C 3 - 16 C^A 3 - E3B3D3 

~2B 2 3 E 3 -4BlA 3 + 4C 3 Bl 
Ai = 2E 3 A 3 ~D 3 B 3 , A 2 = 2 B 2 3 + D 3 B 3 - 2 E 3 A 3 - 8 A 3 C 3 , 
A 3 = 4 BjA 3 + 16 C 3 2 A 3 -4C 3 B 2 3 -D 2 3 C 3 + E 2 3 A 3 + 8 C 3 E 3 A 3 ~ 16 A 2 3 C 3 

-4 C 3 B 3 D 3 + 8 C 3 D 3 A 3 - 4 B 3 E 3 A 3 , 
A 4 = 4 B 2 3 A 3 + 16 C 2 3 A 3 -\ C 3 B 2 3 - D 2 3 C 3 + ElA 3 + 8 C 3 E 3 A 3 -16 A 2 3 C 3 

-4C 3 B 3 D 3 -8C 3 D 3 A 3 + 4B 3 E 3 A 3 . 

(2) 

The expressions g 2 and A5 are polynomials of degree ten in the variables .A3, B 3 , C 3 , D 3 and E 3 
(for instance, g^ contains 147 monomials), g^ is a polynomial of degree 20 in the same variables. 
Their explicit expressions involve too many terms to be reproduced here. 

2 One should note that this is just an auxiliary parametrization and not a uniformization of eq. ([tl|). 
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Expressing the coefficients A3, ■ ■ ■ , E3 in terms of the entries of the i?-matrix, one discovers 
further factorizations: 

g£> = (ci di - c 2 d 2 ) ■ G 2 l > , A 2 = {c t d 1 -c 2 d 2 )-5 2 , 
A 3 = (ci di - c 2 d 2 ) 2 ■ 63, A 4 = (ci di - c 2 d 2 ) 2 ■ 64 

where G 2 is the sum of 1570 monomials of degree ten, S 2 is the sum of 104 monomials of degree 
six, 63 and 64 are the sum of 780 monomials of degree eight, and Ai is the sum of 256 monomials 
of degree eight in the entries of the three-dimensional i?-matrix, i.e. a, ■ ■ ■ , d§. 

Let us note that A 2 = 2 (see equation (|4(])), and is thus related to the projective invariant 
2proj of the conic (|38|). 



5.4 Subcases of 830 

The three-dimensional model B3D was built in such a way that it "projects" down to the two 
dimensional Baxter model, as defined in section ||. It is natural to consider the conditions on B3D 
obtained by writing that the three projections lie on the critical or disorder varieties of the Baxter 
model. 

For example, writing the three disorder conditions as + ds = bs + cb @ of the Baxter model 
for the three projections (i = 1,2,3), yields a codimension-three subvariety of the three-dimensional 
model parametrized as follows: 

a = bi + b 2 + &3 — 2z, c\ = b\ + d\ — z 

c 2 = b 2 + d 2 - z, 03 = 63 + d 3 - z (49) 

On the subvariety (^9|), the discriminant of P/^s') vanishes, A 2 = A4 = 0, the conic ( |38| ) degener- 
ates since 2proj = A 2 = and even more remarkably, Pi{s') gets proportional to (s' — 2) 4 . 

Similarly, the three criticality conditions clb = b& + Cb + ds fbj yield a codimension three 
subvariety: 

2di = a + 61 — 62 — h — 2ci, 2d 2 — a + b 2 - b 3 - bi - 2c 2 , 

2d 3 = a + b3-h-b 2 - 2c 3 (50) 

On the subvariety (|50|), the discriminant of Pi(s') vanishes, A 2 = A 3 = 0, and now Pa(s') gets 
proportional to (s' + 2) 4 . This last codimension-three subvariety is particularly interesting, since 
it is T3D invariant. 



5.5 A solvable case. 

Another interesting model is obtained by setting d\ — d 2 — c?3 = 0. The projections yield six-vertex 
models |26[ . The biquartic equation (|3^) becomes a homogeneous equation of degree 4 and the 
solution is the union of four lines r' — Xr. Remarkably, for this model the left-hand side of ( |35| ) 
factorizes for r' = —r: 

r (ab 2 + 6163 - c 2 + abi + 6 2 6 3 - c\ — 2ac 3 — 26 3 c 3 + 2cic 2 ) 

(ab 2 + 6163 — c 2 + abi + b 2 b 3 — c\ + 2ac 3 + 26 3 c 3 — 2cic 2 ). (51) 

What is more interesting is that we can get the conditions for the existence of solutions to ( fi"l| ) 
in this case. In (|ll|), there are two equations which fix uniquely some scale factors: 

a = /!, apqr — ^ip'qr'. (52) 
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Using these informations, the six others components of (|Tl]) fall into three pairs of equations like: 

ap' = bip + c 3 q + c 2 r, 

a— = 61-+C3-+C2-. (53) 
p' p q r 

In writing these equations, we discarded the trivial solution p = q = r = which always exists in 
this case. Multiplying pairwise these equations, we obtain three linear equations for the variables 
X p = q/r + r/q, X q = r/p + p/r and X r = p/q + q/p. This system of equation can be easily solved. 
These three variables are not independent since they depend only on the ratios of p, q and r. They 
satisfy the relation: 

X p X q X r - {X 2 p + X 2 q + X 2 ) +4 = 0. (54) 

Rewritten in the homogeneous variables a, b\, b 2 , 63, C\, c 2 , C3, equation ( |54| ) is a necessary 
condition for equation ( |Tl"| ) to have non-trivial solutions. This is however not the end of the story 
since the relation pqr = p'q'r' yields another condition on R once we have solved for r/p and q/p. 
Note that the normalization of the variables p, q, r and p' , q' and r' remains free. 

A complete analysis therefore yields the existence of non trivial solutions to ([ll]) when R is on 
some codimension-two subvariety in the parameter space. 



6 Conclusion. 

We have shown how to associate algebraic curves with three-dimensional vertex models. We have 
described a specific model for which the analysis of these curves is handable. We have introduced a 
generalization of the quadratic Frobenius relations (associated to elliptic functions). It corresponds 
to new intertwining relations of products of more than two algebraic curves by i?-matrices living 
on algebraic varieties which are no longer curves. In the example detailed in this paper, one has 
an intertwining of three curves by an i?-matrix living on a highcr-dimcnsional algebraic variety. 
We think that these equations are a key ingredient for the construction of the generalization of the 
Bethe Ansatz in higher dimensions, the quest of solutions of the tetrahedron equations and more 
generally any exact calculation (inversion trick quest of critical manifolds psf ) performed on 
higher dimensional models. 
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